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We present a gauge-invariant effective action for the Abelian Higgs model (scalar electrodynamics) 
with a chemical potential p on a 1+1 dimensional lattice. This formulation provides an expansion 
in the hopping parameter k which we test with Monte Carlo simulations for a broad range of the 
inverse gauge coupling /3 p i and small values of the scalar self-coupling A. In the opposite limit of 
infinitely large A, the partition function can be written as a traced product of local tensors which 
allows us to write exact blocking formulas. Their numerical implementation requires truncations 
but there is no sign problem for arbitrary values of p. We show that the time continuum limit of 
the blocked transfer matrix can be obtained numerically and, in the limit of infinite f3 p i and with 
a spin-1 truncation, the small volume energy spectrum is identical to the low energy spectrum of a 
two-species Bose-Hubbard model in the limit of large onsite repulsion. We extend this procedure for 
finite /3 p i and derive a spin-1 approximation of the Hamiltonian. It involves new terms corresponding 
to transitions among the two species in the Bose-Hubbard model. We propose an optical lattice 
implementation involving a ladder structure. 

PACS numbers: 05.10.Cc,05.50.+q,11.10.Hi,11.15.Ha,64.60.De,75.10.Hk 


I. INTRODUCTION 

The lattice formulation of quantum chromodynamics 
has provided successful treatments of nonperturbative 
problems involving strongly interacting particles in the 
context of particle and nuclear physics. However, nu¬ 
merical computations at finite density and real time have 
remained a major challenge. For this reason, there has 
been a lot of interest in the possibility of building quan¬ 
tum simulators for lattice gauge theory (LGT) using op¬ 
tical lattices (for recent reviews see Refs. p][2j). 

When physical degrees of freedom - individual cold 
atoms or condensates trapped in the optical lattice - play 
the role of the gauge fields, one needs to make sure that 
the physical observables are gauge invariant. One way 
to achieve this goal is to impose Gauss’s law at least in 
some approximation. This question is discussed in the 
context of Abelian gauge theories in Refs. [3H9]. An¬ 
other approach that has been advocated (10, m is to 
work directly with a gauge invariant formulation. 

A simple example is the correspondence between the 
Fermi-Hubbard model and a SU(2) gauge theory with 
one fermion m■ At strong gauge coupling and small 
hopping parameter, the quadratic part of the gauge- 
invariant effective Hamiltonian in one spatial dimension 
reads 

H eff <x 1 + 2 (BjBi +1 + B\ +1 Bi )] , (1) 

i 

where Mi and Bi are the gauge invariant operators cor¬ 
responding to mesons and SU (2)-baryons. A similar 
effective Hamiltonian can be obtained at second order 
in degenerate perturbation theory of the Fermi-Hubbard 
model with strong onsite repulsion. 

Effective actions for mesons and baryons can be con¬ 
structed for SU(N) gauge groups with similar approxi¬ 
mations m ■ The computation of the corrections due to 


the plaquette interactions is rather involved EH- Vari¬ 
ous related techniques have been developed to approach 
chiral symmetry breaking at strong coupling [13 [Hi and 
can accommodate a chemical potential mi. Despite this 
progress, testing the validity of the approximations with 
numerical methods has remained challenging which sug¬ 
gests to consider simpler examples. 

For the Bose-Hubbard model with one species of parti¬ 
cle, there is a remarkable level of quantitative agreement 
OH between state of the art quantum Monte Carlo (MC) 
calculations and their experimental optical lattice imple¬ 
mentations. It would be highly desirable to provide a 
similar proof of principle for a simple LGT model where 
calculations at finite density and real time are possible. 

In this article, we explain how to establish an ap¬ 
proximate quantitative correspondence between the well- 
known Abelian Higgs model (scalar electrodynamics) 
on m with a chemical potential p on a 1+1 dimen¬ 
sional lattice and specific many-body theories that can 
in principle be realized experimentally on optical lat¬ 
tices. The model describes a charged complex scalar 
field and an Abelian gauge field. Its formulation is re¬ 
viewed in Sec. El The model can be studied numeri¬ 
cally by various numerical methods in various regions of 
the parameter space and a first important step is to con¬ 
nect the standard isotropic space-time calculations to the 
Hamiltonian approach obtained by taking the time con¬ 
tinuum limit first and then connecting this Hamiltonian 
with the Hamiltonian of cold atoms trapped in optical 
lattices. Our first goal is to achieve some calibration in 
the gapped, confining phase corresponding to the Mott 
insulator phase of the cold atom system. 

The exact form of the gauge-invariant effective action 
for the mesons is derived in Sec. |III| by integrating over 
the gauge fields. The absence of fermions permits in¬ 
expensive MC simulations at Euclidean time and zero 
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chemical potential. This allows us to test the expansion 
of the effective action in the hopping parameter for a 
broad range of the inverse gauge coupling /3 p i and small 
values of the scalar self-coupling A. Good agreement with 
MC is shown in Sec. ITVl 

In Sec. [V]we consider the opposite limit of arbitrarily 
large A where the amplitude of the scalar held is frozen 
to unity. In this limit, all the remaining variables in¬ 
volved in the basic formulation are compact and we can 
use recently developed Tensor Renormalization Group 
(TRG) methods to write the partition function 

as a traced product of local tensors with discrete indices. 
This reformulation allows us to write exact blocking for¬ 
mulas which can be used for numerical purposes. The 
practical implementation requires truncations but in con¬ 
trast to the MC approach there is no sign problem for 
arbitrary values of f.i . We then write the transfer matrix 
as a product of tensors along a time slice. 

The time continuum limit of the blocked transfer ma¬ 
trix can be obtained numerically. In the limit of infinite 
/3 p i and with a spin-1 truncation (inspired by gauge mag¬ 
net or gauge link constructions fTS, 1 ,[2B]), the small volume 
energy spectrum is identical to the low energy spectrum 
of a two-species Bose-Hubbard model in the limit of large 
onsite repulsion mm- We calculate numerically the 
energy spectrum in this limit in Sec. m 

In Sec. VlTJ we extend this procedure for finite (3 p i and 
derive a spin-1 approximation of the Hamiltonian. It in¬ 
volves new terms corresponding to transitions among the 
two species in the Bose-Hubbard model which we discuss 
in Sec. |VIII| We then propose an optical lattice imple¬ 
mentation involving a ladder structure and discuss fur¬ 
ther plans to obtain a good correspondence between con¬ 
ventional MC calculations and optical lattice measure¬ 
ments. 


II. THE ABELIAN HIGGS MODEL 

In this section, we briefly remind the reader of the 
action for the Abelian Higgs model on a 1+1 space-time 
lattice of size N s x N r and introduce the notations used 
later. We use x, y etc. for space-time vectors, i, j etc. 
for the one-dimensional spatial sites and u= s or f for 
the unit vectors in space and time, respectively. The 
gauge fields U x .o = exp iA x ^ are attached to the links. 
We denote the product of C/’s around a plaquette U p i tX 
where x is the lower left corner of the plaquette in space- 
time coordinates. We use the notation f3 p i = 1/e 2 for 
the inverse gauge coupling and k s ( k t ) for the hopping 
coefficient in the space (time) direction. For the potential 
for the complex scalar field (f> x = \(j) x \ exp (iO x ), we follow 
the convention of Ref. [29]. The action reads: 

S = S g + S h + S x , (2) 

where the gauge part is 

Sg = -p p lY / Re l U pl,x}’ ( 3 ) 

X 


the hopping 

S h = -Kr Y [e^x+^+f + Ul,r <t>* 

X 

~~ K S + / > x+s^x,s^’x (4) 

X 

and the self-interaction 

S\ = \Y (^x -1) 2 + Y^*- ( 5 ) 

X X 

The partition function can then be written as 

Z = J D$D<\>DUe~ s . ( 6 ) 


In the hopping part of the action Sh, we can separate the 
compact and non-compact variables 

Sh — 2k t | <J) X 11 c| / ( cos(Q x -\-f 0 X + A xx +/) 

X 

‘2Hs\4 > x\\4 > x+s\ ^Y COs(^x+s @x + ^Xjs)- (7) 
x 

This equation makes clear that the chemical potential is a 
constant imaginary gauge field in the time direction and 
that the Nambu-Goldstone fields 9 X can be eliminated by 
a gauge transformation 

Ax,i> t A x j) + 9 X , (8) 

which leaves the plaquette terms unchanged. 


III. A GAUGE-INVARIANT EFFECTIVE 
ACTION 


As explained in the introduction, unlike other ap¬ 
proaches H SHUE! [30H32] we will not try to implement 
the gauge field on the optical lattice, but rather try to 
implement a gauge-invariant effective action obtained by 
integrating over the gauge fields. In this section, we 
will show that this effective action is a function of the 
composite, gauge-invariant, meson field which we denote 
M x = <$ x <px in order to emphasize the analogy with Eq. 
([I]). In other words, 



For this purpose, we use the Fourier expansion of the 
Boltzmann weights in terms of the modified Bessel func¬ 
tions I n , for instance, 

exp[2ft T | /+11 ^x+t | cos($a;_j--r 0 X + A x ^ x ^/r)] (10) 

oo 

= y ' -^n(2s+ \4 > x\ I/++TI) GXp[?'n($ x _{_'f Ox + A x x 2/i)], 

n =—oo 
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and similar expressions for the space hopping and the 
plaquette interactions. We can then collect all the expo¬ 
nentials involving a given A x p and perform the integra¬ 
tion over A x ,o- This results in Kronecker deltas relating 
the various Fourier modes. The final result is 


-s. 


Sf 


= £ 


{™n} 


Irnu (Ppl) (4«,s (2ft s |0x| |0x+s |) 


. (11) 


X (2/^r \&x | l^x-t-T I) ^ x P(M^G,f ) 


with specific rules to express the link indices n Xi p in terms 
of the plaquette indices mu that we now proceed to ex¬ 
plain. Given that for x real, I n (x) = /_„( x) there are 
several equivalent ways to label the contributions. We 
use the convention where time is along the vertical axis. 
Starting from the lower left corner of the plaquette and 
moving counterclockwise, the gauge fields in the expo¬ 
nentials come with a plus sign for the first two links and 
a minus sign for the last two links. For the hopping terms, 
the gauge fields always come with a positive sign. This 
results in the rules 


^below 'W'above ( 12 ) 

T^X-f = 'bright ^left (13) 


where the subscripts such as “below” refers to the pla¬ 
quette location with respect to the link. This completely 
fixes the link indices in terms of the plaquette indices 
and it is easy to check that Eqs. ( I2p3 | guarantee that 
the link indices automatically satisfy the current conser¬ 
vation imposed by the integration of the 9 X variables. In 
other words, the mu are the dual variables Illl- 

Equations ( Ppl ) have simple electromagnetic 
analogs. First, n Xi f can be interpreted as a charge 
and m x as an electric field in the spatial direction. 
With this Minkowskian interpretation, Eq. (13) enforces 
Gauss’s law. Second, n xt p can be interpreted as a 
two-dimensional current and m x as a magnetic field 
normal to the two-dimensional plane. In this Euclidean 
interpretation, Eqs. (12 13) express the current as the 
curl of the magnetic held. A right hand rule can be 
obtained for the following index ordering: time, space, 
normal direction. A discrete version of Gauss’s theorem 
guarantees that the sum of the charges on a time slice is 
the last m on the right minus the first to on the left. 

At the lowest order of the strong-coupling expansion 
we have fi p i = 0 and from /„(0) = 0 for n ^ 0, we see that 
all the indices must be zeros. The effect of the plaquette 
can be restored perturbatively. This can be organized in 
an expansion in the hopping parameter. In the isotropic 
case k t = n s = k we obtain: 


S -ff = J2 i-^MxM y + ^K 4 (M x My)A (14) 
(xy) V ' 

- 2k4 T77TA £ M x M y M z M w + 0 (k 6 ) 

l0 ^ l) U^vzw) 



FIG. 1. 1/0 at dpi = 20 for A = 0.05 and A = 0.1 as function of 
k compared with the hopping expansion at 0 p i = oo at 0(k 3 ) 
and 0 (k 5 ). 


IV. MONTE CARLO CALCULATIONS 

Consider the action ([2| for the isotropic case k t = 
k s = n. We start with the /3 P z —> oo limit when all gauge 
variables are frozen to unity, U Xj p = 1. In this limit the 
expectation value of the hopping term 

£* = <Re{4 U x ,p<f> x+0 }) (15) 

can be calculated with the hopping parameter expansion, 
for small k. It has been derived up to 0(k 5 ) in Ref. 291 . 
This result can be generalized to /3 < oo by including the 
appropriate factors of /i(/3)// 0 (/3) for the diagrams that 
involve plaquettes [55] , 

To check the range of validity of the expansion we per¬ 
form Monte Carlo simulations at several values of /3, n 
and A on a 16 2 lattice. To test the /3 —I► oo limit we 
set /3 = 20 and for A = 0.05 and 0.1 scan the range of 
k € [0.05,0.30]. The results for L ^ are shown in Fig. [T] 
The lines represent the expansion at two orders. The ex¬ 
pansion starts to break down around k = 0.15 at 0 (k 3 ) 
and k = 0.2 at 0 (k 5 ). At the present, exploratory, stage 
we use the updating algorithm of Ref. |34| that is appli¬ 
cable at small self-coupling A. Monte Carlo tests of the 
model at large A are left for the future. 

To study the dependence on /3 we focus on the k £ 
[0.08,0.16] range at A = 0.1 and perform calculations at 
Ppl = 20, 2, 0.2 and 0.02. The results for L^ are shown 
in Fig. [2] 

To understand the dependence on /3 p i better we also 
calculate for several values of (3 p i at fixed k = 0.15. 
The results together with the hopping expansion are 
shown in Fig. [3] 


V. THE LARGE A LIMIT 

We now turn to the limit where A becomes arbitrar¬ 
ily large. In this limit, M x is frozen to 1, or in other 
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FIG. 2. L 0 at /3pz = 20, 2, 0.2 and 0.02 for A = 0.1 as func¬ 
tion of k compared with the hopping expansion with included 
dependence on /3 p i up to 0(k 5 ). 
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FIG. 3. at fixed k = 0.15 as function of f3 p i for A = 0.1 
compared with the hopping expansion with included depen¬ 
dence on /3 p i up to 0(k 5 ). 


FIG. 4. The basic B and A tensors (in brown and green, re¬ 
spectively, colors online). The W S 1 are associated with the 
vertical tensors, and the horizontal (spatial) links of the lat¬ 
tice. The A ^ are associated with the horizontal tensors, and 
the vertical (temporal) links of the lattice. 


Following the general principles of the construction [22], 
we attach a B tensor to every plaquette 


£(□) 

1711171217131714 


(19) 


tm D (M, 

o, 


if mi = m 2 = m 3 = 7714 = mu 
otherwise, 


a Al( s ) tensor to the horizontal links 
4( s ) — 1 , 

m a bovembeiow L | mfceZoiu — mabo 


.|(2/e a ), (20) 


and a A^ tensor to the vertical links 


words, the Brout-Englert-Higgs mode becomes infinitely 
massive. We are then left with compact variables of inte¬ 
gration in the original formulation ( 9 X and A X y) and the 
Fourier expansions described before lead to expressions 
of the partition function in terms of discrete sums. 

As explained in Ref. [22], these sums can be formulated 
in a compact way using tensorial notations. In order to 
simplify the equations below we define 

tn(z) = (16) 

Hz) 

These normalized Bessel functions have useful properties. 
For instance t n ( 0) = S n ^. For z non-zero and finite, we 
have 

1 = t 0 (z) > ti(z) > t 2 (z) > ■ ■ ■ > 0. (17) 

In addition, for sufficiently large z, 

t n (z)* 1 -g. (18) 


A ( T ) 

mieftlTlright 


= t\m left -m Tight |( 2 « t ) e< m — 

( 21 ) 


The partition function can now be written as 
Z = {HP p i)H2k s )I 0 {2k t )) v x (22) 


Tr 


TT 

11 ^rriat 


AV BW 

rni e ftm r ight 1711171217131714 


The traces are performed by contracting the vertical and 
horizontal indices as shown in Fig. [4] Note that the ten¬ 
sor associated with a horizontal space link is repre¬ 
sented by a vertical line orthogonal to the link and the 
tensor A^A associated with a vertical time link is repre¬ 
sented by a horizontal line orthogonal to the link. 

The traces can also be expressed in terms of a trans¬ 
fer matrix T which can be constructed in the following 
way. First, we define a matrix B as the product of the 
plaquette tensors on a time slice alternating with the 
link tensor corresponding to the vertical links in between 
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the plaquettes. There are two natural ways to impose 
boundary conditions. The first is to connect the last A 
tensor with the first B tensor (periodic boundary condi¬ 
tions), the second is to impose m = 0 for the first and 
last B tensor (open boundary conditions). However, in 
both cases, the total charge on the interval has to be 
zero. A more general option is to allow arbitrary to’s at 
each end. This would allow us to consider appropriately 
selected charge sectors. In the following we will focus on 
the open boundary conditions and define 

) = (2^r tmi (fipl ) X 

^|mi— yn 2 \— m3| (2^-7-) ■ • ■ 

^ m Ns (fipl)tm Ns (2 k t ). (23) 

Note that with this choice of open boundary conditions, 
the chemical potential has completely disappeared. If 
we had chosen different m’s at the end allowing a total 
charge Q inside the interval, we would have an additional 
factor exp(/xQ). We next define a matrix A as the prod¬ 
uct 


^(mi,m2,...mN„)(rn' 1 , r n' 2 ...m'N ) (^4) 

^|mi —mj| ( 2 K s)t\m 2 -m ' 2 | (2^s) • ■ A\m Ns —m' N | (2^ s ) ■ 

With these notations we can construct a symmetric 
transfer matrix T. Since B is diagonal, real and positive, 
we can define its square root in an obvious way and write 
the transfer matrix as 

T = V^AvS. (25) 

With this definition, the partition function can be written 
as 


Z = (I 0 (P p 1 )I 0 (2k s )I 0 (2 Kt )) v Ti- [T^] . (26) 

Alternatively, we could diagonalize the symmetric matrix 
A and define the (dual) transfer matrix 

f = VamVI. (27) 


The A and B matrices can be constructed by a recur¬ 
sive blocking methods similar to those discussed in Ref. 
[22] . We can construct a new B tensor, called B ', by 
contracting two B tensors on both sides of a A^ tensor. 
This process can be iterated and is illustrated in Fig. [5] 
Using tensorial notation we write 


O/ _ 

M (mi , 1712 ) M' (m^ ,m' 2 ) 


E 


r? . a{ t ) u . 

L 'm 3 m 4 mim 1 ^ L m 4 m5 ? 


(28) 


where the notation M(mi,m 2 ) stands for the product 
state, M = mi®m 2 . One can continue blocking horizon¬ 
tally as above until the desired spatial size is achieved re¬ 
sulting in the matrix B. By their very nature the B' s en¬ 
force that this object is diagonal in the collective product- 
state indices of the upper and lower tensor legs. For prac¬ 
tical reasons, truncation methods need to be introduced 
at the beginning and after each blocking mm- 



FIG. 5. Part of the construction of the blocked B' tensor. 
This shows the contraction of the B and A i - T> tensors. The 
dashed lines are the links of the original lattice. 


m 1 


m 2 


FIG. 6. Graphical representation of the blocking of the A 
tensors. The vertical tensors are the and the dashed 
lines are the links of the original lattice. 


Now consider the A 1 - 8 ) matrices which correspond to 
the horizontal links of the original lattice as in Fig. [6] 
If one takes their outer product and collects the upper 
indices into a single product-state index, and does the 
same for the lower indices, one has another matrix built 
out of the vertical A^ tensors. Once again in 
notation this can be written as 

4 ' 0 ) _ A (s) A (s) 

M(m\ ,7712)-M' (m^ : m' 2 ) mim'j m2m f 2 ’ 

One can continue taking the product of 
until the desired spatial size has been reached 
in the matrix A. 


tensorial 


(29) 

matrices 

resulting 


VI. THE LIMIT dpi 00 

In this section, we discuss the case where both A and 
dpi are infinite. In this limit, the correspondence be¬ 
tween the classical lattice model and a two-species Hub¬ 
bard model has been outlined in Ref. m and we will 
show that the method discussed in Sec. El allows us to 
match the spectra at small volume with very good accu¬ 
racy. This section is a warm-up for the more complicated 
situations discussed in Secs. IVIII and IVIIII 

In the limit where f3 p i becomes arbitrarily large, we 
need to maximize R eU p i iX . This is accomplished when 
A Xy p is a gauge transform of 0. In other words, we can 
set A Xy 0 = 0 and the gauge equivalent configurations 
are taken into account by integrating over the Nambu- 
Goldstone modes 0 X . We recover the familiar 0(2) model. 
As explained in Ref. |36h38| . the time continuum limit 
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can be achieved by taking the limit n T —> oo while keep¬ 
ing the product k t k s constant. This leads to a Hamil¬ 
tonian for quantum rotors located at each spatial site 
and having quantized angular momentum running over 
positive and negative integers. In order to realize this 
Hamiltonian on optical lattices, we considered [27] a spin- 
1 truncation where the angular momentum at each site 
is restricted to the values 0 and ±1. This truncation has 
very small effects on the phase diagram provided that 
the hopping parameter and the chemical potential are 
not too large. 

Following Refs. [271128] with the replacements /?„ —> 
2k 1/ , the Hamiltonian for the spin-1 approximation of the 
0(2) model reads: 

" = (£<<>) : ’ ( 30 > 

i 

-a£ £ ‘<>- lEtWw+Ww) • 

i i 

with the now dimensionful quantities 


U = 


1 - _ M T — 2K S 

X ? h — 5 ^ ; 

Zn T a a a 


(31) 


and a the time lattice spacing. If we use open boundary 
conditions, the sums run over all spatial sites in the first 
two sums and all but the last in the third sum. This 
approximate Hamiltonian can be matched UDdEI with 
the effective Hamiltonian of a Bose-Hubbard model dis¬ 
cussed in sec. VIII that, we expect, can be implemented 
on optical lattices. 

The first step is to understand the matching in the 
simple situation where U is large compared to p and J. 
In this case, U sets the scale of the mass gap and p and 
J introduce small splittings. We now use energy units 
where U= 1. The time continuum limit can be obtained 
numerically as we keep blocking the transfer matrix in 
the space direction as discussed in Sec. [V] The spectra 
in blocks of size 2, 4 and 8 in the spin-1 approximation are 
shown in Fig. [7] for J=0.1, p = 0. A nonzero chemical 
potential would split the charge conjugated states and 
make the graph difficult to read. An example will be 
shown for L = 2 in Sec. ECED 


VII. THE TIME CONTINUUM LIMIT AND 
THE ENERGY SPECTRUM 


In this section, we construct the time continuum limit 


of the transfer matrix T defined by Eq. (251 starting 


now with the more general situation where f3 p i is finite. 
In order to obtain Hamiltonians corresponding to Bose- 
Hubbard models, a guiding strategy is to assume that 
the spatial hopping k s is proportional to some hopping 
energy. In the limit k s = 0, the inside matrix A becomes 
the identity and the transfer matrix T = a/BA-v/B be¬ 
comes B which is diagonal. In this limit, the only way 


Spectra for L=2, 4 and 8; J/L/=0.1; /v=0 



FIG. 7. 0(2) spectra in U units for L= 2, 4, and 8, with 

J=0.1, p = 0. Some higher energy states are not shown on 
the figure. 


to obtain a time continuum limit, in other words to have 
T close to the identity, is to require that both k t and (3 p i 
become large. At leading order in the inverse of these 
large parameters, the eigenvalues of T are 


1 -+ - 
2 p p l 


(32) 


L N. 


) + 


1 

2 K, 


(m 2 + (to 2 - mi ) 2 + ... 

• • • + ( m Ns - m Ns _i) 2 + to 


JV S 


There are two limiting situations: 1 << j3 p i « n T and 
1 « n T « j3 p i. In the first case, 1//3 p i is the largest 
coefficient and the low energy states are those with a 
few m.i nonzero and consequently only a few differences 
of to’s are nonzero and the second term in Eq. (32) is 


a perturbation. In the second case, 1 /k t is the largest 
coefficient and we could guess that to minimize its con¬ 
tribution, we need to take most of the to’s to be equal 
and just create a few charges. For a large volume, this 
can cause the first term of Eq. (32) to dominate, which 


signals an infrared instability. In electrostatic terms, cre¬ 
ating one charge also creates a constant electric field over 
the entire volume. The linear potential is responsible for 
the confinement and trying to treat the interaction per- 
turbatively results in infrared problems. 

For this reason, in the following, we will only consider 
the case 1 << f5 p i « n T and set the scale with the 
(large) gap energy 


U P = 


af3 P i 


(33) 


It is important to d istinguish this scale from the U in¬ 
troduced in Sec. VI Up is associated with the plaquette 
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quantum number m (the energy to create a flux tube 
across one lattice spacing between the opposite charges), 
while U in Sec. VI is associated with the quantum num¬ 


ber n (the energy to create a single charge). In addition, 
we define the (small) energy scales 


Y = 


1 _ Ppi jy p 

2 n T a 2 k t ’ 


which plays the same role as U in Sec. 
hopping parameter 

X = V2(3 p in s Up. 


VI 


(34) 


and the space 


(35) 


We are now in position to derive an expression for the 
Hamiltonian in the spin-1 approximation where the pla- 
quette quantum number m takes values ±1 and 0 only. 
The effect of k s can be studied by linearization. In the 
case of two spins, we find that 


OT/cbc s k =0 = V2 (L x {1) + Lf 2] ) (36) 


We use the notation LX to denote the first generator of 
the spin-1 rotation algebra at the site (1). The notation 
L is used to emphasize that the spin is related to the m 
quantum numbers attached to t he q laauettes in contrast 


to the spin-1 generators L in Sec. VI having a spin related 


to the charges n attached to the time links. The final 
form of the Hamiltonian H for 1 << /3 p i « n T is 


B =^'E( L t»y+ < 37 > 

i 

f<). 

i i 


where )+ ' is a short notat ion to include the single terms 
at the two ends as in Eq. (32), i.e. besides (LX — L^) 2 , 
(L z {2 ) ~ ^( 3 )) 2 ’ ’ (i(jv s -i) - L Us) ) 2 terms this sum 

contains (LX ) 2 and (L Z N )) 2 - 


The process outlined in this paper can in principle be 
generalized for non-Abelian gauge groups in 1+1 and 
higher spatial dimensions. Tensor formulations and ex¬ 
act blocking procedures have already been worked out for 
non-Abelian gauge groups in Ref. [22]. These questions 
are under active consideration. 


VIII. A TWO-SPECIES BOSE-HUBBARD 
MODEL IMPLEMENTATION 


In Refs. [27j|28], we have proposed to match the Hamil¬ 
tonian of the 0(2) model given in Eq. (30) with the two- 
species Bose-Hubbard Hamiltonian on a square optical 
lattice 


H=- y^(t Q a|fflj + t b b\bj + h.c.) - ^(p a+b + A a )nf 
(ij) i,a 

+ TT- !) + W Y 1 n i n i + V » n t n< j ( 38 ) 

i,a i (ij)** 


with a = a, b indicating two different species and with 
nf = a\a,i and n\ = b\bi . In this expression, the chem¬ 
ical potential p a +b is associated with the conservation 
of n a + n b and should not be confused with the chem¬ 
ical potential introduced in the previous section which 
couples to n a — n b and breaks the charge conjugation 
symmetry. In the limit where U a = U b = U and W and 
fJa+b = (3/2)17 are much larger than any other energy 
scale, we have the condition nf + n\ = 2 for the low en¬ 
ergy sector. The three states |2,0), 11,1) and |0, 2) satisfy 
this condition and correspond to the three states of the 
spin-1 projection considered above. 

Using degenerate perturbation theory Pi ns , we 
found [23 [28] that 



+ -jr 7 + L i L t) + (U 0 - IV) 5>r) 2 

° (ij) i 


+ 

+ 




p(n + l)t 2 

b u 0 




(39) 


where p is the number of neighbors and n is the occupa¬ 
tion (p = 2, n = 2 in the case under consideration). L 
is the angular momentum operator in the representation 
n/ 2. The effective hamiltonian (39) without the last term 
has been studied in m and shows a rich phase diagram. 

In order to match with the 0(2) model, we need to tune 
the hopping amplitude as t a = \JV a U /2 and have J = 
VKH, 0 = 2 (U-W), and fx = -(A a -V a ) + (A b -V b ). 


-(A< 

The matching is illustrated in Fig. Nfor J/U = 0.1 and 
large U. In the limit U —» oo, the matching is excellent. 
As we lower XJ, high energy states get a lower energy and 
can mix with the low energy states. These states may 
be related to the modes that become infinitely massive 
when A —> oo. 

The optical lattice implementation is discussed in Ref. 
m for a 8 'Rb and 41 K Bose-Bose mixture where an in¬ 
terspecies Feshbach resonance is accessible. The inter¬ 
species interaction (IV) can therefore be controlled by an 
external magnetic field. The extended interaction, V a , 
is present and small when we consider Wannier gaussian 
wave functions sitting on nearby lattice sites [42] , 

For the Hamiltonian H in Eq. (371 corresponding to 
1 << P p i « k t . we need to introduce a new interaction 
to represent the L x effects that interchange the m = 0 
states with the m = ±1 states. This can be achieved by 
adding the piece 


AH = -^^ 2 ^+^ ■ ( 4 °) 

i 


The matching between the two models can be achieved by 
imposing t = 0, V a = V b = —Y/2 and t ab = X. Energy 

















0(2) and BH Spectra for L=2; JIU=0A; */=0.02 


Abelian-Higgs and BH Spectra for L=2; XIU P =YIU P = 0.1 



0(2) and BH Spectra for L=4; J!U=QA\ {j=0 



Abelian-Higgs 


Abelian-Higgs and BH Spectra for L=4; XIU P =YIU P =0A 


Abelian-Higgs 


FIG. 8. 0(2) with J/U = 0.1 and Bose-Hubbard spectra for FIG. 9. Abelian-Higgs model with X/U = 0.1, Y/U = 0.1 
L = 2 with ft, = 0.02 (top) and L = 4 with jl = 0 (bottom). and the corresponding Bose-Hubbard spectra for L = 2 (top) 

and L = 4 (bottom). 


levels obtained numerically with this method have a good 
matching shown in Fig. [9j Note that if an interspecies 
nearest neighbor interaction V a b is introduced, it can be 
incorporated in 77 e // by making the substitution V a + 
Vb —> V a + Vb — 2V ab in the first term of Eq. (391. 


This is a very different realization than for the 0(2) 
limit. The presence in the Hamiltonian of the additional 
term that interchanges the species index, Eq. (37), rules 
out implementing this Hamiltonian using mixtures of two 
different types of atoms, and it also makes it difficult to 
realize it using two hyperfme states of the same atomic 
species. It could be realized with a single atomic species 
on a ladder structure with a and b corresponding to the 
two legs of the ladder. Ladder systems have been real¬ 
ized experimentally by using lattices of double wells m 
147] . The hopping amplitudes can be tuned such that the 
hopping in the direction along the ladder is negligible, 
but finite along the rungs, thus exchanging a, b species 
index at the same rung. An attractive intraspecies in¬ 
teraction (V a = V b = —y/2) is also needed, favoring 
having two atoms in neighboring sides on the same leg 
of the ladder. This acts as a nearest-neighbor ferromag¬ 
netic coupling in the effective one-dimensional spin chain 
Hamiltonian. For the experimental implementation, an 
attractive nearest neighbor interaction can be obtained 
by using cold dipolar atoms or molecules, with dipole mo¬ 


ments aligned along the ladder and with inter-rung dis¬ 
tance such that the rapidly decaying dipole-dipole inter¬ 
action between next-nearest-neighbors can be neglected. 
With this alignment of the moments, the inter-species in¬ 
teraction in the same rung (W) is repulsive. Experiments 
with ultra cold dipolar quantum gases have been per¬ 
formed with chromium [ 45] , erbium SSI ISO], and dyspro¬ 
sium [5lfl53| . which have magnetic moments 6/is, 7/is 
and 10/is, respectively, and with polar molecules, such 
as 40 K 87 Rb and 23 Na 40 K [HH5S]- 


IX. CONCLUSIONS 


With the remarkable experimental progress in the field 
of Atomic, Molecular and Optical physics, particularly 
the loading of trapped ultra-cold atoms onto optical lat¬ 
tices, and the unprecedented levels of control and tunabil- 
ity of their interactions, a wide range of quantum many- 
body lattice Hamiltonians can be realized and studied. 

We have proposed a two-species Bose-Hubbard model 
that may be used as a quantum simulator for lattice 
gauge theories. The gauge invariance is built-in and 
thus does not need to be achieved via fine-tuning, and 
the correspondence between the proposed Bose-Hubbard 
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Hamiltonian and the Abelian Higgs model can be checked 
quantitatively. 

On the experimental side, various multi-species Bose- 
Hubbard systems can be created with current cold atom 
technology. Various possible realizations for the two 
species that are needed can be explored. On the the¬ 
ory side, the Bose-Hubbard Hamiltonian is amenable 
to Quantum Monte Carlo calculations and the Abelian 
Higgs model can be treated with the TRG (and the time- 
continuum limit obtained) as well as Monte Carlo calcu¬ 
lations. The full spectra can be calculated, at least for 
small systems. The TRG formulation can very naturally 
provide time-dependent correlation functions. Therefore, 
beside the ultimate goal of quantum simulating lattice 
gauge theories using cold atoms on optical lattices, the 


correspondence we report here may also provide a way of 
studying the dynamics of the cold atom system such as 
response to sudden quenches. 
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